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Abstract 

The momentum distribution of itinerant electrons in the one-dimensional 
Falicov-Kimball model is calculated for various ground-state phases. In par- 
ticular, we examine the periodic phases with period two, three and four (that 
are ground-states for all Coulomb interactions) as well as the phase separated 
states (that are ground states for small Coulomb interactions). For all periodic 
phases examined the momentum distribution is a smooth function of k with 
no sign of any discontinuity or singular behavior at the Fermi surface k = kp- 
An unusual behavior of (a local maximum) is found at k = 2>kF for elec- 
tron concentrations outside half-filling. For the phase separated ground states 
the momentum distribution exhibits discontinuity aX k = k^ < kp. This 
behavior is interpreted in terms of a Fermi liquid. 

PACS nrs.: 71.27.+a, 71.28.+d, 71.30.+h 
Keywords: Falicov-Kimball model, momentum distribution, Fermi liquid, Luttinger 
liquid 
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1 Introduction 



Since its introduction in 1969 the Falicov-Kimball model PP has become an important 
standard model for a description of correlated fermions on the lattice. The model 
describes a two-band system of localized / electrons and itinerant d electrons with 
the short-ranged f-d Coulomb interaction U. The Hamiltonian is 

H = Y: U.dfd, + UY: fthdtdi + Ef Y.fth (1) 

ij i i 

where f^, ft are the creation and annihilation operators for an electron in the local- 
ized state at lattice site i with binding energy Ef and df, di are the creation and 
annihilation operators for an electron in the conduction band. The conduction band 
is generated by the hopping matrix elements tij, which describe intersite transitions 
between the sites i and j (as usually, we assume that = —t ii i and j are nearest 
neighbors and tij = otherwise, and all energies are measured in units of t). 

The model has been used in the literature to study a great variety of many-body 
effects in metals, of which valence and metal-insulator transitions, charge- density 
waves and electronic ferroelectricity are the most common examples |21 E] • It has 
been applied to a variety of lattices, one [H E], two [HI 13 El , three 0, and infinite 
dimensional fU] , and occasionally to small clusters [TTl IT^ IT^ . Exact results are 
available in very few instances fOlElElEl and general theorems have been proved 
for special cases [7j. In spite of the existence of an analytic solution in o? = oo dimen- 
sion ^1 and an impressive research activity in the past, the properties of this 
seemingly simple model are far from being understood. Indeed, while the ground- 
state properties of the /-electron subsystem has been satisfactory explained, only a 
few exact results are known concerning the ground-state correlations of the itinerant 
electrons ^Hj- Even, some important quantities such as a momentum distribution of 
itinerant electrons has not been explored yet. The first attempt to describe ground- 
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state correlations of itinerant electrons has been performed recently by Jedrzejewski 
et al. • They calculated a number of one and two-point correlation functions in or- 
der to describe short and long-range correlations between itinerant electrons, as well 
as between itinerant and localized electrons for various ground states. Calculations 
have been done analytically and by means of well-controlled numerical procedures. 
In this paper we use the same method to calculate a momentum distribution of itin- 
erant electrons for physically the most interesting cases. In particular, we examine 
the periodic phases with period two, three and four as well as the phase separated 
states that are ground states for small (large) /-electron concentrations. 

2 The method 

Before showing results for the momentum distribution of itinerant electrons in the 
one-dimensional Falicov-Kimball model let us briefly summarize the main steps of 
the calculational method Since in the spinless version of the Falicov-Kimball 
model without hybridization the /-electron occupation number f^fi of each site i 
commutes with the Hamiltonian (1), the /-electron occupation number is a good 
quantum number, taking only two values: = 1 or 0, according to whether or not 
the site i is occupied by the localized / electron and the Hamiltonian (1) can be 
written as 

^ = E tijdfdj + f/ E ^^dtd^ + EfNf, (2) 

ij i 

where Nj = J2i Wi denotes the number of / electrons. 

Thus for a given /-electron configuration w = {wi, W2 ■ ■ ■ wl} defined on a one-di- 
mensional lattice (of L sites) with periodic boundary conditions, the Hamiltonian (2) 
is the second- quantized version of the single-particle Hamiltonian h{w). If we denote 
by {|«)}i=i,...,L the orthogonal basis of one-electron states, such that df creates an 
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electron in the state then the matrix elements of h{w) in the basis {|«)}i=i,...,L 
are defined by 

H{w) = jz mw)\j)dtd, + EfNf, (3) 
where the non-vanishing matrix elements are given by 

{i\h{w)\i) = Uwi, {i\h{w)\j) = -1 if 3 = i±l. (4) 

Let {|z^)}j/=i/i,...,,y^ be the orthonormal basis built out of the eigenstates of h{w) to 
the eigenvalues X^, such that X^, < Xpi if < v' . Then, the unitary matrix U, with 
the following matrix elements Uiy-. 

Uiu = (5) 

diagonalizes the matrix of h{w), 

$]W+(j|/iH|z)W,., = AA.'. (6) 

hi 

Moreover, the set of operators {6^, b^}, z/ = i/i, . . . , z/^,, defined by 

L 

1=1 

satisfies the canonical anticommutation relations, and 

H{w) = mw)\j)dtd, + EfNf = XXK + EfNf. (8) 

i,j=l y 

Since, 

d, = J2iMiy)K, (9) 
we can express the momentum distribution function of itinerant electrons 

n, = -^Ee^'^'"'^(4rf^) (10) 
in terms of the site-components l\i\v) of the eigenvectors \v): 

n^ = \Y.^'''^'~'^ E(Hj)(^k), (11) 

j,Z v<vp 



where up stands for the label of eigenvectors |z/), such that there is exactly A^^ 
eigenvectors with u < up. For some low-period ion configurations the eigenproblem 
can be solved exactly, while for any other ground-state configurations of interest 
the numerical exact-diagonalization procedure can be used. In the next section we 
calculate the momentum distribution of itinerant electrons for three periodic phases 
with the smallest periods as well as for the phase separated ground states. Such a 
selection of phases is not accidental. In our previous papers fSlElEOl we have shown 
that just these phases occupy the largest regions in the {Ef — U) ground-state phase 
diagram of the Falicov-Kimball model. Indeed, for large U there is only one nontrivial 
ground-state phase (besides the empty or fully occupied phase), and namely the two- 
period phase w^^^ = {10 . . . 10}, and for intermediate values of U there are only two 
other relevant phases and namely the three-period phase w^'^^ = {110. . . 110} and 
the four-period phase w^^^ = {1110. . .1110}. These phases are ground states also 
for small values of U, but in this region also the phase separated configurations, 
as well as the periodic phases with larger periods can be the ground states of the 
Falicov-Kimball model. The periodic phases are insulating, and the phase separated 
configurations are metallic and thus one can expect fully different behavior of the 
momentum distribution in these ground states. 

3 Results and discussion 

As was described above one has to know the site-components of the eigen- 
vectors |z/) to calculate the momentum distribution of itinerant electrons. In order 
to calculate {ilu) for w^^\ w^"^^ and w^^^ we have generalized the procedure used 
by Jedrzejewski et al. ^Hj for analytical calculations of correlation functions in the 
two-period phase. 
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For w^'^\ n = 1, 2, 3 the nonvanishing matrix elements of h(w^^^) are given by 



+ 1) = -1, {i\h{w^''^)\t) = Uwt^. 



(12) 



One can easily verify that the matrix of h{w^^'^) can be rewritten to a block- diagonal 
form by reordering the original basis k = 2'n-l/L, / = 0, 1, . . . , L — 1. The new 

basis for w^'^\ n = 1, 2, 3 that reduces h{w^^'^) to a block-diagonal form is given by: 



{\k),\k + 7r)}, n = l, 

{|A;),|A; + 27r/3),|fc + 47r/3)}, n = 2, 

{|A;), \k + 7r/2), |A; + vr), |A; + 37r/2)} , n = 3, 



(13) 



where k = 2nl/L, I = 0,1, . . . , L/{n + 1) - 1. 

The diagonal blocks are n + + 1 matrices h{w^"'^) that in the corresponding 

L J k 

basis given by preceding equation have the form 

L Jfc -a Ek+n + a 



U 



(14) 



/i(i/;(2)y 



Ek + 2a —a —a 
-a £fc+27r/3 + 2a -a 
-a -a efc+47r/3 + 2a 



U 



a 



(15) 



efc + 3a —a —a 

—a eA:+7r/2 + 3a —a 

—a —a £^fc+7r + 3a 

—a —a 



-a 
-a 
-a 



-a Ek+3^/2 + 3a 
where Ek = —2 cos k. 

Let, X^^^k) be eigenvalues of h{w^"'^), i.e., 

M^^"^)] J^)?^ = ^'';\k)\k)'-;\ /i = 1, 2, . . . n + 1. 
Then the corresponding eigenvectors lA;)^"-* can be written as: 



U 



a 



(16) 



(17) 



l^^(i) ^ y(i)|A; + vr), 

1^)^ = + y^k + 27r/3) + |A; + 47r/3), (18) 

= x^^^k) + |fc + 7r/2) + z(^)\k + vr) + + 37r/2), 
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where explicit expressions for coefficients x^^\y'^^\ z'^^\v^^\ n = 1,2,3 are given in 
Appendix. 

Since we are interested in the half-filled band case Nd + Nf = L (which is the 
point of the special interest for valence and metal-insulator transitions caused by 
promotion of electrons from localized / orbitals to the conduction band states) we can 
restrict our considerations to the lowest A^^^ eigenvalues Xi^\k) and the corresponding 
eigenvectors \k)i^\ and set 

Then, the site-components (m|/c) for considered phases w^"'\n = 1,2,3 read 
{m\ky^^ = ('x(^)e^*^'" + ?/(i)e*('=+^>'") , 

(^|A;)(2) = M2)e^'=™ + 2/(2)e*('=+M3)™. ^ ^(2)g^(fc+4V3)m^ ^ ^20) 

(m\k)^^^ — _J_ frf>{S)^ikTn _^ y{S)^i{k+7v/2)m _^ ^(3) ^i{k-\-7v)m _^ y{3) ^i{k-\-37v/2)m\ 

Substituting these expressions into (fTT|) and doing some tedious algebra one ob- 
tains the final expression for the momentum distribution of itinerant electrons 

n« = |x»r, ^ = 1,2,3. (21) 

Using this expression (and corresponding expressions from Appendix) one can 
plot the momentum distribution of itinerant electrons in the particular phase as a 
function of k for different values of the Coulomb interaction IJ . For the two-period 
phase w^^^ the results obtained are displayed in Fig. 1. It is seen that for all examined 
V the momentum distribution is a smooth function of k. There is no sign of any 
discontinuity at the Fermi surface k = kp = t\ j1 (the no n- Fermi liquid behavior) as 
well as no sign of singular behavior aX k = kp (the non-Luttinger liquid behavior). 
This can be verified analytically since the expression ()21|) (for i = 1) can be rewritten 
as 
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Obviously, there is no discontinuity as well as singular behavior at kp for finite 
U. Near kp, n^f}^ behaves like 



what corresponds exactly to the behavior of the Tomonoga-Luttinger fermions cou- 
pled by 2fci? potential [23 • This potential gives rise to a gap in the energy band 
spectrum of the Tomonoga-Luttinger fermions, and the appearance of this gap re- 
sults in the smooth behavior of the momentum distribution. We believe that the 
same mechanism (the existence of a gap in the charge excitation spectrum at half- 
filling) is responsible for a smooth behavior of the momentum distribution of the 
Falicov-Kimball model. The same behavior at half-filling exhibits also the Hubbard 
model that has a gap in the charge excitation spectrum at n-f = = ^, but not 
for any other density. Since the Falicov-Kimball model can be considered as an ap- 
proximation to the full Hubbard model in which one part of electrons (say with spin 
down is immobile), it could be interesting to compare our solution ()22|1 with known 
results for the momentum distribution of the Hubbard model. 

The analytical results for the momentum distribution nj^ in the Hubbard model 
are known in the strong coupling limit where the first two terms of the perturbation 
expansion read |23 




(23) 




(24) 



Comparing this result with the strong coupling expansion of ()22|) 




(25) 



one finds that the Fahcov-Kimball model (in spite of the above mentioned simplifi- 
cation) contains still much of physics of the full Hubbard model at half-filling. This 
result is not surprising since for both models there is the antiferromagnetic long-range 
order in the ground state for = n/ = n-f = = |. Since this ground state per- 
sists for all Coulomb interactions U > (for both the Hubbard and Falicov-Kimball 
model) one could expect a good accordance of results also for intermediate and weak 
interactions. For arbitrary U there is only an approximate (an antiferromagnetic 
Hartree-Fock) solution for nj! in the full Hubbard model ? that has precisely the 
some form as our solution for the Falicov-Kimball model. 

Outside the half-filling, the momentum distribution behaves, however fully dif- 
ferently in the Falicov-Kimball and Hubbard model. While in the Hubbard 
model [23] exhibits the power-low singularity at k = kp for rif = 7^ | (even 
for U 00), Uk in the Falicov-Kimball model remains still a continuous function of 

k for all finite Coulomb interactions. The situation for 71^ = ^ (the phase w^'^^) is 

(2) 

displayed in Fig. 2. A continuous character of is obvious for all finite Coulomb 
interactions U, and increasing U only smears . In accordance with the Hubbard 
model an unusual behavior of the momentum distribution is observed near = Skp 
and U small. However, while the momentum distribution has a weak singularity at 
k = Skp in the Hubbard model it has a local maximum in the Falicov-Kimball 

model at this point. The similar behavior is observed also for nd = ^ (see Fig. 3). 

(s) 

Again there is a continuous change of at /c = kp, with no sign of any discou- 
rse 

tinuity or singular behavior for finite U. Also an unusual behavior of n\. (a local 
maximum) is observed at k = 3k p. 

To reveal the origin of this unusual behavior we have performed the lowest order 
perturbation calculations of ?7,^*\ i=l,2,3 in terms of U. A straightforward application 
of the perturbation procedure jSHI yields the following expression for the momentum 
distribution of itinerant electrons in the Falicov-Kimball model (up to the second 
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order) : 



n 



,T/(0 12 

(i) _^2j2 \ k,k'\ (0)^ ioT \k\>kF, (/ = 1,2,3), (26) 



k' {^k - £k') 

where 



<i' = fEe^^'=-'^'^^^^f (27) 

and nf^ is the Fermi distribution function of noninteracting electrons. 

For the periodic phases w'^'-' (1=1,2,3) the matrix elements V^''l, can be directly 
calculated and the off-diagonal elements, that enter to Eq. (|^. are given by 



y(0 ^ I -ITT' k' = k±2nkF, n = l,2.../, ^^g^ 



(0 



I 0, otherwise. 
Substituting this expression into Eq. one obtains the final expression for 

^(0 = I ^ (29) 

for kpil + 2{n - 1)) < k < kpil + 2n). 

Analysing this expression one can find that n^^^ changes its behavior at points 
k = 3kF,5kp,7kF, . . . and this change can produce the local maxima at some points. 
Indeed, we have found that the second order perturbation theory can describe the 
local maximum at A; = 3A;i? for Ud = ^- This is illustrated in Fig. 4 where the 
momentum distribution (exact as well as perturbation results) is displayed for the 
periodic phases with the period three, four and six. It is clear that this unusual 
behavior of the momentum distribution at /c = Skp, bkp, Ikp ... is caused by the 
periodic arrangement of the localized / electrons, which results in a very simple 
form of off-diagonal matrix elements V^^\,. It is interesting that a similar unusual 
behavior of the momentum distribution at k = Skp, 5kp, 7kp . . . has been observed 
(predicted) also in the Hubbard model, however in this case it was assigned to a pair 
(multipairs) of electron-hole excitations 
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Finally, let us briefly discuss the case of phase separation. It is well-known 
that the ground states of the Falicov-Kimball model for sufficiently small Coulomb 
interactions {U < 1) and sufficiently small {rid < 1/4) or large (3/4 < rid < 1) 
d-electron concentrations are the phase separated configurations, i.e., configurations 
in which one-half of lattice is fully empty or fully occupied by /-electrons. Such 
configurations are metallic and we expect fully different behavior of the momentum 
distribution of itinerant electrons in this phase. Unfortunately, the phase separated 
configurations are not periodic and thus we cannot proceed in the analytic calcula- 
tions as in the preceding cases, however, the numerical procedure is still possible. To 
calculate numerically in the phase separated region {U and small) we need to 
know exactly the ground-state configuration w^^ for selected U and n^. In general, 
it is very difficult to find the ground-state configuration for arbitrary U and Ud from 
the phase separated region, however, it is possible to solve this task for some spe- 
cial values of U and n^. For example, exhaustive small-cluster exact- diagonalizat ion 
studies that we have performed in our previous paper ^3] for the Falicov-Kimball 
model in the weak-coupling limit showed that the ground-state configuration of the 
model for U = 0.6 and Ud = 1/8 is of the type wP' = {111100 ... 111100111 ... 111}. 
The momentum distribution n^^ for this configuration (calculated numerically us- 
ing the procedure described in the preceding section) is displayed in Fig. 5a. As 
we conjectured the momentum distribution of itinerant electrons behaves fully dif- 
ferently in the metallic phase. It seems that there is a discontinuity (or a singular 
behavior) in n^'^ at some critical value oi k = < kp- To determine exactly 
which type of behavior realizes near k^ we have calculated the finite-size discontinu- 
ity A = nP/^ 2n — n^J 2-n as a function of 1/ L. The results obtained are plotted in 
Fig. 5b and they clearly show that there is a discontinuity at /c = k^. This shows on 
the Fermi liquid behaviour of itinerant electrons (or some part of them) in the phase 
separated (metallic) state. To verify this conjecture we have calculated separately 
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contributions to n^^ from different parts of a lattice. Since the configuration can 
be formally considered as a mixture of two phases w"' = {111100 ... 111100} and 
w'' = {111. ..111} it is natural to divide the lattice into two parts a and b and to 
rewrite the summations in Eq. (Illj) as follows 

< = ^ E e''^'-'^ E (Hj) {iW) = E = < + < + <^ + < 

j,l V<Vp j,l 

= E E E E ^.v- (30) 

j&a,l^a j&a,l^b j€b,l&a j&b,l&b 

The momentum dependence of single contributions nl^n^^^n^^ and n\ is plotted 
in Fig. 6a for two values of L. It is seen that with increasing L the contribution n^^ 
+ n^" goes to zero and thus in the thermodynamic limit L — oo only nl and n\. 
remain finite. The first contribution ra^ exhibits the same behavior as the momentum 
distribution in the periodic (insulating) phases w'^^\ w^"^^ and w^^\ There is no sign 
of any discontinuity or singular behavior at k = kp. This result is expected since also 
the phase w"" is periodic. In contrary to this case the second contribution n\ exhibits 
an unexpected behavior of the Fermi liquid type with discontinuity at some critical 
momentum /c = /cq = O.lvr that does not coincides however with the Fermi surface 
momentum kp = ^vr = 0.1257r. Such a behavior as well as the meaning of k^ can 
be easily explained within the formalism described above (a mixture of two phases). 
Let A^^ be the number of itinerant electrons in the phase w^. It can be calculated 
directly from the expression for (Eq. (fTT|l ) putting j = I and taking the sum only 
over the lattice sites from w'^, i.e, 

iv.' = E E(Hj)(j»- (31) 

jgfe V<Vp 

In Fig. 6b we plotted the quantity kp = -jf'^, that represents the Fermi surface 
of itinerant electrons in w^, as a function of 1/L (L^ = 5L/8 denotes the size of w'^ 
phase). It is seen that kp goes to ko for L — oo what provides a clear physical 
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interpretation for ko. According these results the itinerant electrons in the phase 
behaves like the Fermi liquid with a discontinuity at kp = jfir = kg. 

In summary, the momentum distribution ra^ of itinerant electrons in the one- 
dimensional Falicov-Kimball model has been calculated for various ground-state 
phases. In particular, we have examined the periodic phases with period two, three 
and four (that are ground-states for all Coulomb interactions) as well as the phase 
separated states (that are ground states for small Coulomb interactions). We have 
found that for all periodic phases examined the momentum distribution is a smooth 
function of k with no sign of any discontinuity or singular behavior at the Fermi 
surface k = kp- An unusual behavior of Uk (a local maximum) is found at k = Skp 
for electron concentrations outside half-filling. For the phase separated ground states 
the momentum distribution Uk exhibits discontinuity aX k = ko < kp. This behavior 
has been interpreted in terms of a Fermi liquid. 
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4 Appendix 

Here we give the explicit expressions for coefficients x^^\y^^\ z^^\v'^^'^ from (fTHj) and 
the explicit expressions for eigenvalues A^"')(fc) of h{w^"'^), (n = 1,2,3): 



(1) ^ I q(l) y(l) = 



,(1) 



1-6' 



(1) 



(32) 
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(2) _ iz£i^ (2) (2) _ lz££^ (2) (2) _ (2) 



(3) _ lz4^o(3) y(3) _ 1-4'^ q(3) ^(3) _ 1-4'^ g(3) ^(3) _ (3) 



where 



bl\k) =a(^\k + 2n/3), c^^\k)=a(^\k + An/3), 



ap{k) = -^e,-3 + ^X(^Kk), 

hf{k)=af{k + 7^/2), c(f)(A:) = a(f)(fc + vr), {k) = a^^) {k + 3ix /2), 



and 



A5'^(A:) = f + vW+4, 

\f\k) = -2rcos(|) + f 
Af (A;) = 2rcos(f + |) + f 

(A;) = 2rcos(f - |) + 

q = cos{3k) + (f )3, r = ^(f )2 + 1, = arccos (^) , 

Ap)(A;) = -^p - ^(48a - 6g - 288cg^i - 2Aa'^q-^ + 72v^6j9-i)i/2 + |?7, 
A(^)(fc) = -^p + ^(48a - 6g - 288cg-i - 24a2g-i + 72^/Qbp-^y/^ + |f/, 



Xf{k) = '§p- ^(48a -6q- 288cq-^ - 24a\~^ - 72^/Qbp^^f''^ + |[/, 
Ai^^fc) = ^p+ ^(48a - 6g - 288cg-i - 2Aa^q~^ - 72^hp~^fl^ + f f/, 



a = 3f/V8 + 4, b = Uy8, c = f/V4 - 3f/7256 - 2 cos(4A:) + 2, 
p = (4a + g + 48cg~^ + 4a2g~^)^/2, 
g = (288ac + 1086^ - 8a^ + 12s)l/^ 

s = (-768c3 + 384c2a2 - 48ca^ + 432ac62 + 816^ - Ub'^a^y/'^. 
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Figure Captions 



Fig. 1. The momentum distribution of the one dimensional Falicov-Kimball model 
calculated for = | and several different values of U. 

Fig. 2. (a) The momentum distribution of the one dimensional Falicov-Kimball 
model calculated for = | and several different values of U. (b) The weak-coupling 
results for the momentum distribution on the enlarged scale. 

Fig. 3. (a) The momentum distribution Uk of the one dimensional Falicov-Kimball 
model calculated for rid = \ and several different values of U. (b) The weak-coupling 
results for the momentum distribution on the enlarged scale. 

Fig. 4. (a) The lowest-order perturbation results for the momentum distribution 
of the Falicov-Kimball model calculated for three different (i-electron concentrations, 
(b) A comparison of exact and perturbation results obtained for U = 0.1 and two 
different (i-electron concentrations. 

Fig. 5. (a) The momentum distribution n^** of the one dimensional Falicov- 
Kimball model calculated numerically for = |, t/ = 0.6 and L = 3200 sites, 
(b) The finite-size discontinuity A = n^J'.^Tv, — n^i!'_2n as a function of 1/L. 

Fig. 6. (a) Contributions to the momentum distribution n^^ of the one dimen- 
sional Falicov-Kimball model at = | from different parts of lattice {U = 0.6 and 
L = 3200). (b) k'^p = ^TT as a function of 1/L. 
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